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Abstract 

We review the present status of rare B decays, focusing on inclusive decay 
modes and their role in our search for new physics. We also briefly discuss 
direct CP violation in rare B decays and the rare kaon decays -—>■ tt^uu 
and Kl n^uu, which offer complementary opportunities for precision flavour 
physics. 
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1 Introduction 



Flavour physics deals with that part of the standard model (SM) which dis- 
tinguishes between the three generations of fundamental fermions. It is still a 
mystery why there are exactly three generations. Also the origin of the fermion 
masses and their mixing is unknown; in particular, the SM does not explain the 
hierarchical pattern of these parameters. Flavour physics can be regarded as the 
least tested part of the SM. This is reflected in the rather large error bars of 
several flavour parameters such as the mixing parameters at the 20% level 
which has to be compared with errors smaller than 1% in high energy electroweak 
precision experiments. 

However, the experimental situation concerning flavour physics is drastically 
changing. There are several B physics experiments successfully running at the 
moment and, in the upcoming years, new facilities will start to explore B physics 
with increasing sensitivity and within various different experimental settings: 
Apart from the CLEO experiment (Cornell, USA), located at the Cornell Electron- 
Positron Storage Ring (CESR) 0, two B factories, operating at the T(4S') reso- 
nance in an asymmetric mode (fig.|l]), have started successfully: the BaBar experi- 
ment at SLAC (Stanford, USA) and the BELLE experiment at KEK (Tsukuba, 
Japan) . Besides the successfully running hadronic B physics program at FER- 
MILAB (Batavia, USA) ^ there are independent B physics experiments planned 
at the hadronic colliders: the LB.C-B experiment at CERN in Geneva and the 
i?TeV experiment at FERMILAB ||^ . The main motivation for a B physics pro- 
gram at hadron colliders is the huge b quark production cross section with respect 
to the one at e~^e~ machines. 

While the time of the electroweak precision physics focusing on the gauge 
sector of the SM draws to a close with the completion of the LEP experiments at 
CERN and the SLC experiment in Stanford, the era of precision flavour physics 
focusing on the scalar sector of the SM has just begun with the start of the B 
factories. 

The B system represents an ideal framework for the study of flavour physics. 
Since the b quark mass is much larger than the typical scale of the strong interac- 
tion, long-distance strong interactions are generally less important and are under 
better control than in kaon physics thanks to the heavy mass expansion. Thus, 
for example the CP violation in the B system will yield an important independent 
test of the SM description of CP violation (see |§]). B meson decays also allow 
for a rich CKM phenomenology and a stringent test of the unitarity constraints. 

The so-called rare decays are of particular interest. These processes represent 
flavour changing neutral currents (FCNC) and occur in the SM only at the quan- 
tum level. The inclusive rare decay modes are theoretically clean observables 
because no specific model is needed to describe the hadronic final states. Their 
role is twofold: on the one hand they are relevant to the determination of CKM 
matrix elements. On the other hand they are particularly sensitive to new physics 
beyond the SM, since additional contributions to the decay rate, in which SM 
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particles are replaced by new particles such as the supersymmetric charginos or 
gluinos, are not suppressed by additional factors a;/(47r) relative to the SM con- 
tribution. This makes it possible to observe new physics indirectly - a strategy 
complementary to the direct production of new (supersymmetric) particles. The 
latter production is reserved for the planned hadronic machines such as the LHC 
at CERN, while the indirect search of the B factories already implies significant 
restrictions for the parameter space of supersymmetric models and will thus lead 
to important clues for the direct search of supersymmetric particles. It is even 
possible that these rare processes lead to the first evidence of new physics by a 
significant deviation from the SM prediction, for example in the observables con- 
cerning direct CP violation, although it will then be difficult to identify in this 
way the new structures in detail. But also in the long run, after new physics has 
already been discovered, these decays will play an important role in analyzing in 
greater detail the underlying new dynamics. 

Although the general focus within flavour physics is at present on B systems, 
kaon physics offers interesting complementary opportunities in the new physics 
search such as the rare decays n^uu and Kl —>■ n'^uu. They are specifi- 

cally interesting in view of the current experiments at the Brookhaven laboratory 
(USA) and suggested experiments at FERMILAB (USA) and at KEK (Japan). 

The paper is organized as follows: in Section 2 we briefiy discuss the role of 
the strong interaction within fiavour physics. In Section 3 the status of rare B 
decays within the SM is reviewed. In Section 4 we explore the implications of 
these decays for our search of physics beyond the SM. In Section 5 we discuss 
direct CP violation and in Section 6 the complementary role of rare kaon decays 
within precision fiavour physics. In Section 7 we present our summary. 

2 Strong interaction in B decays 

Flavour physics is governed by the interplay of strong and weak interactions. 
One of the main difficulties in examining the observables in fiavour physics is 
the infiuence of the strong interaction. As is well known, for matrix elements 
dominated by long-distance strong interactions there is no adequate quantitative 
solution available in quantum field theory. The resulting hadronic uncertainties 
restrict the opportunities in fiavour physics significantly. The present discussion 
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Figure 2: Penguin decays of B mesons. 



on the new g — 2 muon data also reflects this issue. While the hadronic self- 
energy contribution to the g — 2 observable can be determined by experimental 
data, the well-known light-by-light contribution can only be modelled at present 
(see for example |I0[1). 

However, there are several fundamental tools available, which are directly 
based on QCD. High hopes for precise QCD predictions are placed on lattice 
gauge theoretical calculations. While there are competitive predictions from lat- 
tice gauge theory for form factors of semileptonic decays, pure hadronic decays 
are less accessible to these methods (fill). With the help of the so-called QCD 
sum rules, a consistency test between hadron physics and perturbative QCD, it 
becomes possible to connect hadronic and fundamental QCD parameters directly. 
Theoretical predictions via QCD sum rules, however, always have relatively large 



uncertainties [12|. Another approach is the method of factorization |13|. This 



method has recently been systematized for nonleptonic decays in the heavy quark 



limit [0]. However, within this approach a quantitative method to estimate the 
l/rub corrections to this limit is missing. The latter contributions can be specifi- 



cally large if they are chirally enhanced []T5[. Further fundamental methods whose 



applications and precision are also somewhat restricted are chiral perturbation 
theory |TB| and heavy quark effective theory [jlT . 



In view of this, the goal must be to minimize theoretical uncertainties with 
the help of an optimized combination of different fundamental methods solely 
based on QCD. This can only be done for a selected number of observables in 
flavour physics. However, there are also observables, dominated by purely per- 
turbative contributions, which will make precision flavour physics possible in the 
near future. Among them inclusive rare B decays (see fig. H |T^) play the most 
important role. 

Inclusive decay modes are theoretically clean and represent a theoretical lab- 
oratory of perturbative QCD. In particular, the decay width r{B ^ ^sl) is well 
approximated by the partonic decay rate r(6 57), which can be analysed in 
renormalization group improved perturbation theory: 

r(5 ^ x,7) = r(6 ^ 57) + A"°'^f^''*- (2.1) 

Nonperturbative effects, A"""^^'^*', play a subdominant role and are under control 



3 



w 




Figure 3: QCD corrections to the decay B —>■ Xs^y. 



thanks to the heavy mass expansion. 

Thus, in general, inclusive decay modes should be preferred to exclusive ones 
from the theoretical point of view. The inclusive modes B — ^ Xs(^ii)1 and B 
Xs{d)l'^l~ can be measured by the electron-positron experiments {B factories, 
CLEO) with their kinematic constraints and their low background, while thy are 
more difficult to measure at hadronic machines. Exclusive decay modes, however, 
are more accessible to experiments, in particular at hadronic machines. But 
in contrast to the inclusive modes, they have in general large nonperturbative 
QCD contributions. Exclusive decays such as Bd^s ~^ 1^^ l^ i Bd K*'j and 
Bd K*^~^fi~ are distinguished observables at the LHC--B experiment. 

Within inclusive B decay modes, short-distance QCD effects lead to a tremen- 
dous rate enhancement. These effects are induced by hard gluon exchange be- 
tween the quark lines of the one-loop electroweak diagrams (fig. |^). 

The QCD radiative corrections bring in large logarithms of the form a"(mf,) 
log™' (mfe/M), where M = rrit or M = my/ and m < n (with n = 0, 1, 2, ...). This 
is a natural feature in any process where two different mass scales are present. In 
order to get a reasonable result at all, one has to resum at least the leading-log 
(LL) series 



with the help of renormalization group techniques. Working to next-to-leading- 
log (NLL) precision means that one is also resumming all the terms of the form 



A suitable framework to achieve the necessary resummations of the large logs 
is an effective low-energy theory with five quarks, obtained by integrating out 
the heavy particles, which, in the SM, are the top quark and the W boson. The 
standard method of the operator product expansion allows for a separation of an 
amplitude of a weak meson decay process into two distinct parts, the long-distance 
contributions contained in the operator matrix elements and the short-distance 
physics described by the so-called Wilson coefficients (see fig. ^). In the case of B 
decays, the W boson and the top quark with mass bigger than the factorization 
scale are integrated out, that is removed from the theory as dynamical variables. 
The effective hamiltonian can be written 



<(mb) log"(m,/M), (LL) 



(2.2) 



asinib) Kinib) log" (m^/M) ) , (NLL). 



(2.3) 



eff - 



4Gf 

71" 



(2.4) 
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Figure 4: Operator product expansion. 



where Oi{fi) are the relevant operators and Ci{^, Mheavy) are the corresponding 
Wilson coefficients. As the heavy fields are integrated out, the complete top 
and W mass dependence is contained in the Wilson coefficients. Working out 
a convenient set of quantities, both in the effective (low-energy) theory and in 
the full (standard model) theory, and requiring equality (matching) up to terms 
suppressed by higher powers of mw or rrit, these coefficients can be determined. 

Within this framework QCD corrections for the decay rates are twofold: the 
ingredients are the order corrections to the matrix elements of the various 
operators and the order corrections to the Wilson coefficients, of course both 
at the low-energy scale /i ~ rrif,. Only the sum of the two contributions is renor- 
malization scheme and scale independent; in fact, from the //-independence of the 
effective Hamiltonian, one can derive a renormalization group equation (RGE) 
for the Wilson coefficients Cj(/i): 

f^-^Ci{fx) = jjiCjifi) , (2.5) 

where the matrix 7 is the anomalous dimension matrix of the operators Oi. Then 
there are the following three principal steps leading to the leading-log (next-to- 
leading-log) result within the effective field theory approach: 

• Step 1: One has to match the full SM theory with the effective theory 
at the scale = fiw: where fiw denotes a scale of order rriw or m^. At 
this scale, the matrix elements of the operators in the effective theory lead 
to the same logarithms as the full theory calculation. Consequently, the 
Wilson coefficients Ci{fiw) only pick up small QCD corrections, which can 
be calculated in fixed-order perturbation theory. In the LL (NLL) program, 
the matching has to be worked out at the 0(a°) (0(q;J)) level. 

• Step 2: Then one performs the evolution of these Wilson coefficients 
from fi = nw down to /i = /ivy, where Hb is of the order of mj,. As the 
matrix elements of the operators evaluated at the low scale /ib are free of 
large logarithms, the latter are contained in resummed form in the Wilson 
coefficients. For a LL (NLL) calculation, this RGE step has to be done 
using the anomalous dimension matrix up to order a\ (a^)- 

• Step 3: To LL (NLL) precision, the corrections to the matrix elements 
of the operators (s7|Ci(/i)|6) at the scale fi = fib have to be calculated to 
order (aj) precision. 
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Finally, we stress that the step from the leading (LL) to the next-to-leading (NLL) 
order within the framework of the renormalization group improved perturbation 
theory is not only a quantitative one increasing the precision of the theoretical 
prediction, but also a qualitative one, which tests the validity of the perturbative 
approach in the given problem. 



3 Inclusive decay modes 

3.1 Experimental data on B ^ Xs^ 

Among inclusive rare B decays, the B Xs'y mode is the most prominent 
because it is the only one that is already measured: in 1993, the first evidence 
for a penguin- induced B meson decay was found by the CLEO collaboration. At 
CESR, they measured the exclusive electromagnetic penguin process B K*'y. 
The inclusive analogue B Xg'-f was also found by the CLEO collaboration 
through the measurement of its characteristic photon energy spectrum in 1994 
(see [^]). As this process is dominated by the two-body decay 6 — >• 57, its 
photon energy spectrum is expected to be a smeared delta function centred at 
^ mb/2, where the smearing is due to perturbative gluon bremsstrahlung 
and to the nonperturbative Fermi motion of the b quark within the B meson. 
Only the high part of the photon energy spectrum is sensitive to the rare decay 
B — >■ Xg^. Some lower cutoff in the photon energy has to be imposed in order to 
exclude the background, mainly from the nonleptonic charged current processes 
b cqq' + 7 or 6 ^ uqq' + 7, which have a typical bremsstrahlung spectrum 
that is maximal at -E^ = and falls off for larger values of E^. Therefore only 
the "kinematic" branching ratio for B — > X^'y in the range between E^ = 2.2 
GeV and the kinematic endpoint at E^ = 2.7 GeV could be measured directly. 
To obtain from this measurement the total branching ratio, one has to know the 
fraction R of the B — > Xgj events with E^ > 2.2 GeV. This was done in where 
the Fermi motion of the b quark in the B meson was taken into account by using 
the phenomenological model of Altarelli et al. (ACCMM model) |^ . Using this 
theoreticalinput regarding the photon energy spectrum, the value R = 0.87±0.06 
was used by the CLEO collaboration, leading to the CLEO branching ratio ||2l|] 



BiB ^ X,7) = (2.32 ± 0.57,tat ± 0.35,,,,) x 10"^ (3.1) 

The first error is statistical and the second is systematic (including model 
dependence). This measurement was based on a sample of 2.2 x 10^ BB events. 
In 1999, CLEO has presented an improved measurement which is based 



on 53% more data (3.3 x 10^ events). They also used the slightly wider E^ window 
starting at 2.1 GeV. The relative error drops almost by a factor of y/S: 

B{B ^ X,7) = (3.15 ± 0.35,t,t ± 0.32,^,, ± 0.26w) x 10^^ (3.2) 

The errors represent statistics, systematics, and the model dependence, respec- 
tively. 
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There are also data at the Z° peak from the LEP experiments. The ALEPH 



collaboration p4[ has measured the inclusive branching ratio 



BiHb ^ X,7) = (3.11 ± OMstat ± 0.72,^,) x 10"^ (3.3) 

It should be noted that the branching ratio in ( |3.3| ) involves a weighted average 
of the B mesons and A;, baryons produced in decays (hence the symbol H^) 
different from the corresponding one given by CLEO, which has been measured 
at the T{4S) resonance. High luminosity is more difficult to obtain at higher 
e~^e~ collision energies. Thus, BB samples obtained by the LEP experiments are 
smaller than the one at CESR. The rate measured by ALEPH is consistent with 
the CLEO measurement. 

Recently, CLEO presented a refined preliminary analysis (with an lower pho- 
ton energy cut > 2.0 GeV) H: 



BiB ^ X,7) = (2.85 ± 0.35stat ± 0.22,^,) x 10"^ (3.4) 



and also BELLE has presented preliminary data ||2^ of competitive experimental 
accuracy: 

BiB ^ X,7) = (3.37 ± 0.53stat ± 0.42,,, ± 0.54^) x 10'^ (3.5) 

More accurate data can be expected in the near future. With the expected high 
luminosity of the i?-factories, an experimental accuracy below 10% in the inclusive 
B — > mode appears to be possible. 

The uncertainty regarding the fraction R of the B —>■ Xgj events with > 
2.2 GeV quoted in the experimental measurement, also cited as model depen- 
dence, should be regarded as a purely theoretical uncertainty. As mentioned 
above, the fraction R was calculated in [l^ using the phenomenological model by 



Altarelli et al., where the Fermi motion of the b quark in the B meson is charac- 
terized by two parameters, the average Fermi momentum pf of the b quark and 
the mass rUg of the spectator quark. The error on the fraction R is essentially 
obtained by varying the model parameters pp and nig in the range for which the 
ACCMM model correctly describes the energy spectrum of the charged lepton in 
the semileptonic decays B — > Xcii' and B — > Xuiu, measured by CLEO and AR- 
GUS. In [|19| a first comparison between the calculated photon energy spectrum 
and the one measured by the CLEO collaboration was presented. The (normal- 
ized) measured photon energy spectrum and the theoretical one are in agreement 
for those values of pp and rrig, that correctly describe the inclusive semileptonic 
CLEO data B Xcii' and B — >■ Xuit^', at present, the data from the radiative 
decays is, however, not precise enough to further constrain the values of pp and 
rrig. The best fit between the theoretical and measured photon energy spectrum 
is obtained for pp = 450 MeV and rrig = 0. One should mention that the analysis 



| T9[1 of the photon energy spectrum, in particular the calculation of the fraction R 
in the ACCMM model used by CLEO, does not include the full NLL information, 
which becomes available in the meantime. 
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Different components of the photon spectrum in the B decay, 



Besides this phenomenological model by Altarelh et al., more fundamental 
theoretical methods are available today to implement the bound state effects, 
namely by making use of operator product expansion techniques in the frame- 
work of heavy quark effective theory (HQET). A new analysis along these lines 
was presented in p6|. Unfortunately, the operator product expansion breaks down 
near the endpoint of the photon energy spectrum; therefore, an infinite number 
of leading-twist corrections has to be resummed into a nonperturbative universal 
shape function, which determines the light-cone momentum distribution of the 



6-quark in the B meson The physical decay distributions are then obtained 
from a convolution of parton model spectra with this shape function. At present 
this function cannot be calculated, but there is at least some information on the 
moments of the shape function, which are related to the forward matrix elements 
of local operators. Ansatze for the shape function, constrained by the latter infor- 
mation, are used. In contrast to the older analysis based on the ACCMM model, 
the new analysis of Kagan and Neubert |^ includes the full NLL information. 
Their fraction R = O.TStoji (for the energy cut > 2.2 GeV) is smaller than 
the factor used by CLEO. 

An important observation is that the shape of the photon spectrum is not 
sensitive to physics beyond the SM. As can be seen in fig. |^, all different contri- 
butions to the spectrum (corresponding to the interference terms of the various 
operators involved) have a very similar shape besides the small 8-8 contribution. 
This implies that we do not have to assume the correctness of the SM in the 
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experimental analysis and, thus, a precise measurement of the photon spectrum 
can be used to determine the parameters of the shape function. 

Clearly, a lower experimental cut decreases the sensitivity to the parameters 
of the shape function (or, more generally, the model dependence). With respect 
to this, the ideal energy cut would be 1.6 GeV. However, in this better 
understanding of the ip background would be mandatory. The intermediate ip 
background, namely B ipXg followed hj ip —>■ X'7 is more than 4 x 10~^ in the 
'total' branching ratio. With the present energy cut of 2.1 GeV this contribution 



is suppressed and estimated to be less than 5% . 

Another future aim should be to determine the shape function (and analo- 
gously the parameter of the ACCMM model) by using the high-precision mea- 
surements of the photon energy spectrum. 



3.2 NLL QCD calculations 



As mentioned above, the inclusive decay B —>■ Xg'j is a laboratory for perturbative 
QCD. Nonperturbative effects (see section |3.3| ) play a subdominant role and are 
well under control thanks to the heavy quark expansion. The short-distance QCD 
corrections enhance the partonic decay rate r(6 — > 57) by more than a factor of 
2. The corresponding large logarithms of the form a"(mb) log"^ (mf,/M) , where 
M = rrit or M = rriw and m < n (with n = 0, 1, 2, ...), have to be summed with 
the help of the renormalization group improved perturbation theory as presented 
in section H. 

The effective Hamiltonian relevant to B ^ in the SM reads 



H,ff{B ^ X,7) 



"71" 



A,EQ(/i)0,(/x) 



(3.6) 



where Oi(/i) are the relevant operators, Ci(/i) are the corresponding Wilson co- 
efficients, which contain the complete top- and W- mass dependence (see fig. 
^), and Xt = VtbV^l with Vij, the CKM matrix elements. The CKM dependence 
globally factorizes, because we work in the approximation A„ = (in the case of 
B —>■ Xsj). One neglects the operators with dimension > 6. which are suppressed 
by higher powers of 1/mw. 

Using the equations of motion for the operators, one arrives at the following 
basis of dimension-6 operators: 



Oi = 






0, = 


{cLal^bia) ( 


SlI31,,ClI3) , 


03 = 






04 = 




{ulpIplUlo) + 


0, = 






Oe = 







+ ... + (bLi3l,,h(} 

La) + ■■■ + {bLlSlt^bla 
+ ...+ (hRp-f^JjRp 
i) + ... + (bRl3^f,bRa 
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Figure 6: SM Hamiltonian in the case of 5 — ^ Xs'y. 



Os = igs/lQn^)s^a'^'' im,ifi)Pn + msifi)PL)T^f,bp G^,. (3.7) 

In the dipole type operators O7 and Os, e and F^u {gs and G^,^) denote the 
electromagnetic (strong) couphng constant and field strength tensor, respectively; 
(a = 1,8) denote SU{3) colour generators. 

The error of the leading logarithmic (LL) result was dominated by a large 
renormalization scale dependence at the ±25% level, which alreday indicated the 
importance of the NLL series. By convention, the dependence on the renormal- 
ization scale fib is obtained by the variation mh/2 < fif, < Sm^. The former 
measurement of the CLEO collaboration (see (|3.1| )) overlaps with the estimates 
based on LL calculations, and the experimental and theoretical errors are compa- 
rable. In view of the expected increase in the experimental precision in the near 
future, it became clear that a systematic inclusion of the NLL corrections was 
becoming necessary. Moreover, such a NLL program is also important in order 
to ensure the validity of renormalization group improved perturbation theory in 
this specific phenomenolgical application. 

This ambitious NLL enterprise was completed some years ago. This was a joint 
effort of many different groups (lT|],[§g, [§3], g),g). The theoretical 
error of the previous LL result was substantially reduced to ±10%, and the central 
value of the partonic decay rate increased by about 20%. 

All three steps to NLL precision listed below (|2.5| ) involve rather difficult 
calculations. 



The most difficult part in Step 1 is the two- loop (or order a^) matching 
of the dipole operators Oj and Os- It involves two- loop diagrams both in 
the full and in the effective theory. It was first worked out by Adel and 
Yao As this is a crucial step in the NLL program, Greub and Hurth 
confirmed their findings in a detailed re-calculation using a different method 
[0 . Two further recalculations of this result [^, |36[ were presented in the 
meanwhile, confirming the original results in ||3^. In order to match the 
dimension-6 operators and Og, it is sufficient to extract the terms of 
order (M = m]^, rrit) from the SM matrix elements for b ^ s'y and 

b sg. Terms suppressed by additional powers of rrib/M correspond to 
higher-dimensional operators in the effective theory. In p4[ the finite parts 
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of the two-loop diagrams in the SM were calculated by means of the well- 
known method of aymptotic mass expansions, which naturally leads to a 
systematic expansion of Feynman diagrams in inverse powers of M. 

The order a1 anomalous dimension matrix (Step 2) has been worked out 
by Chetyrkin, Misiak and Miinz ||33|. In particular, the calculation of the 
elements 7^7 and 7^8 (i = 1, 6) in the 0{a1) anomalous dimension matrix 
involves a huge number of three-loop diagrams from which the pole parts 
(in the d — A expansion) have to be extracted. This extraction was simplified 
by a clever decomposition of the scalar propagator. Moreover, the number 
of necessary evanescent operators was reduced by a new choice of a basis of 
dimension-6 operators. Using the matching result (Step 1), these authors 
obtained the NLL correction to the Wilson coefficient C-j{^b). Numerically, 
the LL and the NLL value for C-j{^h) turn out to be rather similar; the 
NLL corrections to the Wilson coefficient C-j{^b) lead to a change of the 
B Xs7 decay rate that does not exceed 6% [0^ . 

It should be stressed that the result of Step 2, in particular the entries 7^7 
and 7j8 (i = 1, 6) of the anomalous dimension matrix to NLL precision, is 
the only part of the complete NLL enterprise which has not been confirmed 
by an independent group. 

Step 3 basically consists of bremsstrahlung corrections and virtual correc- 
tions. While the bremsstrahlung corrections were worked out some time 



ago by Ali and Greub [IS] and were confirmed and extended by Pott |30 
a complete analysis of the virtual two-loop corrections (up to the contribu- 
tions of the four-quark operators with very small coefficients) was presented 
by Greub, Hurth and Wyler [^. The latter calculation involves two-loop 
diagrams, where the full charm dependence has to be taken into account. 
By using Mellin-Barnes techniques in the Feynman parameter integrals, the 
result of these two-loop diagrams was obtained in the form 

(2 \ " 2 
^ log'"^, (3.8) 

where the quantities Cq and Cnm are independent of rric,. The convergence of 
the Mellin-Barnes series was proved; the practical convergence of the series 
(|3.8|) was also checked explicitly. Moreover, a finite result is obtained in the 
limit rric — >■ 0, as there is no naked logarithm of ml /ml. This observation is 
of some importance in the b ^ d'y process, where the w-quark propagation 
in the loop is not CKM-suppressed (see below). The main result of Step 3 
consists in a drastic reduction of the renormalization scale uncertainty from 
about ±25% to about ±6%. The central value was shifted by about 20%. 



In these results are presented also in the 't Hooft-Veltman scheme, 
which may be regarded as a first step towards a cross-check of the complete 
NLL calculation prediction in a different renormalization scheme. 
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(a) 




Figure 7: (a) Typical diagrams (finite parts) contributing to the matrix element 
of the operator O2 at the NLL level (Step3); (b) typical diagram (infinite part) 
contributing to the NLL anomalous dimension matrix (Step2); typical diagram 
(finite part) contributing in the NLL matching calculation shown in fig.^ (Stepl). 



Quite recently, the results of the matrix elements based on the operator O2 
were confirmed by an independent group with the help of the method 
of asymptotic expansions. 



It is clear that many parts of the NLL calculations at the partonic level in 
the case of 6 ^ S7 can be straightforwardly taken over to the cases b — > d'j, 
b — ^ s gluon and b sl^l^ . In the latter case, however, many modifications 
are necessary; in particular the operator basis gets enlarged as will be discussed 
below. 

Combining the NLL calculations of the three steps, the first complete the- 
oretical prediction to NLL precision for the branching ratio of -B — >^ was 
presented in (see also [^): 



i3(fi ^ X,7) = (3.28 ± 0.33) X 10'^ (3.9) 

The theoretical error has two dominant sources. The dependence, which is 
now reduced to about 6%. The other main uncertainty of 5% stems from the 
rric/mb dependence. This first theoretical NLL prediction already included the 
nonperturbative correction scaling with which are rather small (at the 1% 



level) (see section p^) . Surprisingly, these first NLL predictions ([^,||3^) are 
almost identical to the current prediction quoted in ( p.l7|) , in spite of so many 
important additional refinements such as the electroweak two-loop corrections 
and the nonperturbative corrections which will be discussed below. 



3.3 Nonperturbative contributions 

Within the framework of the heavy mass expansion, nonperturbative corrections 
to the branching ratio of decay B Xg'y can be singled out. These contributions 
also apply to the case of the decay B Xd'j and, with some modifications, to 
the case of the decay B — > XJ~^l~ . 
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Figure 8: a)Feynman diagram from which the operator O arises, b) Relevant 
cut-diagram for the (O2, 07)-interference. 



If one neglects perturbative QCD corrections and assumes that the decay 
B —>■ Xsj is due to the operator O7 only, the calculation of the differential decay 
rate basically amounts to working out the imaginary part of the forward scattering 
amplitude T(g): 

T{q) =i J d^x{B\TOi{x)O7{0)\B) exp{iqx) . (3.10) 

Using the operator product expansion for TO^ix) ^7(0) and heavy quark effective 
theory methods, the decay width r{B Xg'j) reads [^, ^ (modulo higher terms 
in the l/m^ expansion): 



1 , 9 



^Zd = - ^^2 , (3.11) 



where Ai and A2 are the parameters for kinetic energy and the chromomagnetic 
energy. Using Ai = -0.5 GeV^ and A2 = 0.12 GeV^ one gets d^d ^ -4%. 
The B — > Xgj decay width is usually normalized by the semileptonic one. The 
semileptonic decay width gets 1/m^ corrections, which are negative; thus, the 
nonperturbative corrections scaling with 1/m^ tend to cancel in the branching 
ratio B{B Xg'-f), and only about 1% remains. 

Voloshin |^ considered the nonperturbative effects when including also the 



operator 02- This effect is generated by the diagram in Fig. ^a (and by the one, 
not shown, where the gluon and the photon are interchanged); g is a. soft gluon 
interacting with the charm quarks in the loop. Up to a characteristic Lorentz 
structure, this loop is given by the integral 



1 rl-x xy 



As the gluon is soft, i.e. k'^,kgk^ ~ A'^'^^ mi,/2 <^ ml, the integral can be 
expanded in kg. The (formally) leading operator, denoted by O, is 

d = ^V,,V:,C,^^^sj,{l-j,)g.G,,be^^^''d'F,^ . (3.13) 
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Then working out the cut diagram shown in Fig. one obtains the nonper- 

to the decay width, which is due to the (O2, C7) 



turbative contribution T 



{0,07) 

B^Xs-y 



interference. Normahzing this contribution by the LL partonic width, one obtains 



n{0,07) 
-pLL 



1 C2 A2 

'9 CVm2 



^ +0.03 



(3.14) 



As the expansion parameter is rrihls^QCD/m'i ~ 0.6 (rather than Ag^^/m^), it 
is not a priori clear whether formally higher order terms in the mc expansion are 
numerically suppressed. More detailed investigations ||4^, 0] have shown that 
higher order terms are indeed suppressed, because the corresponding expansion 
coefficients are small. 

The analogous 1/m^ effect has been found independently in the exclusive 
mode B K*'y in ref. [0. Numerically, the effect there is also at the few per 
cent level. Moreover, the analysis of the l/fn^ effects was extended to the decay 
B XJ+l- 



m 



44, 46 



As was recently emphasized by Misiak H^, an analogous systematic analysis 



of terms like F 



(02,02) 



at first order in is still missing. Rigorous techniques such 



as operator product expansions do not seem to be applicable in this case. 



3.4 Theoretical status 
3.4.1 B Xsi 

The theoretical prediction for the partonic b —>■ s'-f decay rate is usually normal- 
ized by the semileptonic decay rate in order to get rid of uncertainties related with 
the CKM matrix elements and the fifth power of the b quark mass. Moreover, an 
explicit lower cut on the photon energy in the bremstrahlung correction is often 
made: 

Rquark[0) = ^^TT 7^—^ 3.15 

F[6 ^ Xceue] 

where the subscript S means that only photons with energy > {1 — 6)E!^°'^ = 
(1 — 6)^ are counted. The ratio Rquark is divergent in the limit 5—^1 due to the 
unphysical soft photon divergence in the subprocess b s'ygluon. In this limit 
only the sum of F[6 s'-f], F[6 — > sgluon] and F[6 s'ygluon] is a reasonable 
physical quantity, in which all divergences cancel out. In it was shown that 
the theoretical result is rather sensitive to the unphysical soft-photon divergence; 
the choice 6 = 0.90 was suggested as the optimized definition of the total decay 
rate. In the analysis presented in |]^ the limit 5 ^ 1 is taken and the singularities 
are removed by adding the virtual photon corrections to 6 ^ sgluon. 

It is suggestive to give up the concept of a total decay rate of 6 57 and 
compare theory and experiment using the same energy cut as CLEO {E^ > 
2.1 GeV). Then also the theoretical uncertainty regarding the photon energy 
spectrum mentioned above would occur naturally in the theoretical prediction. 
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In the meanwhile detailed studies of the electroweak corrections were per- 
formed. In 1^ part of the electroweak two-loop contributions, namely contri- 



butions from fermion loops in gauge boson propagators (7 and W) and from 
short-distance photonic loop corrections, were calculated. Moreover, it was found 
that the on-shell value of the fine structure constant l/oem = 137 is more ap- 
propriate for real photon emission than the value l/oem = (130.3 ± 2.3) used 



in previous analyses. The QED loop calculations in |Q confirmed this expecta- 
tion. This change in leads to a reduction of 5% in Rquark- In [0] the QED 
analysis made in was improved by resumming the contributions of order 
a log(/ib/M)(Q;s log(/i{,/M)" to all orders (while in |^ only the n = contribu- 



a 



tion was included). This resummation decreases the QED corrections. In |^ 
complete analysis of the heavy top and the heavy Higgs corrections in the limit 
my/ was made. This analysis was recently refined in [Q. A 2% reduction of 
the branching ratio of -B ^ X^'j due to purely electroweak corrections is found. 

Using the measured semileptonic branching ratio Bf^p , the branching ratio 
B{B — i> Xg'j) is given by 

B{B ^ X,7) = Rg^^rk X BlipXl + Anonpert), (3.16) 

where the nonperturbative corrections scaling with 1/ml and 1/wi^, summed in 
Anonpert, havc a numerical effect of +1% ^ and +3% |^, respectively, on 
the branching ratio only. 

For a comparison with the ALEPH measurement (|3.3|) the measured semilep- 
tonic branching ratio B{Hi, Xc^^^u) should be used consistently. This leads to 
a larger theoretical prediction for the LEP experiments. 

Including only the resummed QED corrections and the nonperturabtive cor- 
rections discussed in section |3.3| , using the on-shell value of aem and working 
with the convention 5 ^ 1 in Rquark, one ends up with the following theoretical 



prediction for the B Xg-y branching ratio |51|: 



B{B X,7) = (3.32 ± 0.14 ± 0.26) x 10-^ (3.17) 

where the first error represents the uncertainty regarding the scale dependences, 
while the second error is the uncertainty due to the input parameters. In the 
second error the uncertainty due to the parameter mc/mi, is dominant. 

Quite recently, quark mass effects within the decay B —>■ X^j were further 
analysed ||55|, in particular the definitions of the quark masses rric and ruh in the 



matrix element {O2) = Xs'y\{sc)v-A{cb)v-A\b) ■ Since the charm quark in the 
matrix element {O2) are dominantly off-shell (see fig. |^) the authors argue that 
the running charm mass should be chosen instead of the pole mass. . The latter 
choice was used in all previous analyses |^ ^ ^ 



m 



P°i7mf^ ^ m^^ifi)/mr'', fxe[me,m,]. (3.18) 



Since the matrix element starts at NLL order and, thus, the renormalization 
scheme for rric and rrih is an NNLL issue, one should regard this choice as an 
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educated guess of the NNLL corrections. However, this new choice is guided 
by the experience gained from many higher-order calculations in perturbation 
theory. Numerically, the shift from rri^°^'^ / m\"^^ = 0.29 ±0.02 to rri^^ [ii) / m\°^^ = 
0.22±0.04 is rather important and leads to a +11% shift of the central value of the 



B Xg'f branching ratio. The authors of ||5^ quote a weighted experimental 
world average using the preliminary data from CLEO and BELLE, ( ^.41 ) and 
(PD, and the published ALEPH data (p|): 

B{B X,7) = (2.96 ± 0.35) x 10"^ (3.19) 

With their new choice of the charm mass and with 6 = 0.9, their theoretical 
prediction for the 'total' branching ratio is 

B{B X,7) = (3.73 ± 0.30) x 10-^ (3.20) 

which means that the difference between the theoretical and the experimental 
value is consistent with zero at the level of 1.6(j (if one assumes that a statistical 
interpretation of this difference is really possible). Because the choice of the 
renormalization scheme for rric and rrib is a NNLL effect, one could argue for a 
larger theoretical uncertainty in m^^{fi)/ml"^'^ which includes also the value of 
cffi^poie^ A more conservative choice would then be m^^(/i)/m^°'^ = 0.22 ± 0.07 
which would reduce the significance of the perceived discrepancy. 

Instead of making a theoretical prediction for the branching ratio B{B 
Xs7), one can use the experimental data and theory in order to directly determine 
the combination | Vj6V^*|/| V^fe| of the CKM matrix elements; in turn, one can 
determine \ Vts\ by making use of the relatively well known CKM matrix elements 
Vcb and Vth- An update of the analysis in |56] was presented in ||5l|]. Using the 



CLEO data (|3), the ALEPH data d^), and the theoretical prediction ( pTTD , 

\v:sVt,\ 



one finds 51 



\v:sV4 



\Vcb\ 

The average of the two measurements yields 

\v:sVtb\ 



0.95 ± 0.08e:.p. ± 0.05t/,. CLEO 
0.91 ± 0.15exp. ± 0.04t,,. ALEPH. 



IK 



cb\ 



0.93 ± 0.09 ± 0.03 = 0.93 ±0.10 (3.21) 



where in the very last step the theoretical and experimental errors were added 
in quadrature. Using |V^b| = 0.99 ± 0.15 from the CDF measurement and \ Vcb\ = 



0.0393 ± 0.0028 extracted from semileptonic B decays, one obtains 



I V"t, I = 0.037 ±0.007, (3.22) 

where all the errors were added in quadrature. This is probably the most direct 
determination of this CKM matrix element. With an improved measurement of 



16 



B{B Xs'j) and Vtt, one expects to reduce the present error on \ Vts\ by a factor 
of 2 or even more. 

Finally, some remarks on the decay mode b — s gluon are in order. The 
effective Hamiltonian is the same as in the b s'-f case. By replacing the photon 
by the gluon, the NLL QCD calculation of b ^ sj can also be used. However, in 
the calculation of the matrix element of the operator O2, further diagrams with 
the nonabelian three-gluon coupling have to be calculated [^. Numerically, one 
obtains B{b —>■ s gluon) = (5.0 ± 1.0) x 10~^, which is more than a factor of 



2 larger than the former LL result B{b s gluon) = (2.2 ± 0.8) x 10" p9 
The mode b s gluon represents one component to the inclusive charmless 
hadronic decays, B — >• Xnocharm, where Xnocharm denotes any hadronic charmless 
final state. A measurement of the corresponding branching ratio would allow 
the extraction of the ratio iKife/Kbl, which is poorly known at present [Q. At 
the quark level, there are decay modes with three-body final states, b q'q'q 
{q' = u,d,s; q = d, s) and the modes b qg, with two-body final-state topology. 
The component b ^ sg of the charmless hadronic decays is expected to manifest 
itself in kaons with high momenta (of order mf,/2), owing to its two-body nature 



54 



3.4.2 B ^ 

With respect to new physics, also the B X^ff decay is of specific interest, 
because its CKM suppression by the factor |VJdP/|V^sP in the SM may not be 
true in extended models. Moreover, a future measurement of the B —>■ X^j 
decay rate will help to drastically reduce the currently allowed region of the 
CKM-Wolfenstein parameters p and 77. 

Most of the theoretical improvements carried out in the context of the decay 



B — * Xs7 (see sections 2]2| and p.3|) can straightforwardly be adapted for the decay 
B Xd'j. As for the former decay, the NLL-improved and power-corrected decay 
rate for B — »• X^j has much reduced theoretical uncertainty, which would allow 
a more precise extraction of the CKM parameters from the measured branching 
ratio. 

The perturbative QCD corrections in the decay B X^i'j can be treated in 
complete analogy to the ones in the decay B —>■ Xg'j. The effective Hamiltonian 
is the same in the processes b ^ sj and 6 — >• ^7 up to the obvious replacement 
of the s-quark field by the (i-quark field. However, as for 6 — >■ 0^7 is not small 
with respect to Xt and Ac, one also has to encounter the operators proportional 
to Xu- The matching conditions Ci^myy) and the solutions of the RG equations, 
yielding Ci{fib), coincide with those needed for the process B Xg^y. The power 
corrections in 1/m^ and l/ni^ (besides the CKM factors) are also the same for 
the two modes. 

The long-distance contributions from the intermediate w-quark in the penguin 
loops, however, are different. These are suppressed in the B Xsj mode by the 
unfavourable CKM matrix elements. In i? ^ -^cz7, there is no CKM suppres- 
sion and one has to include the long-distance intermediate w-quark contributions. 
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which can only be modelled at present. However, these contributions are esti- 



mated to be rather small ^Tj. Moreover, it must be stressed that there is no 
spurious enhancement of the form log(m„/ /if,) in the perturbative contribution to 
the matrix elements (Xrf7|0j„|-B) {i = 1,2) as shown by the explicit calculation 
in [31 1 and also discussed in [58|. In other words, the limit m„ — can be taken. 

The predictions for the B — > X^'j decay given in |^ show that for /x^ = 2.5 
GeV (and the central values of the input parameters) the difference between the 
LL and NLL results is ~ 10%, increasing the branching ratio in the NLL case. 
For a fixed value of the CKM-Wolfenstein parameters p and i], the theoretical 
uncertainty of the branching ratio is: 

ABiB ^ Xd^)/B{B ^ Xrf7) = ±(6 - 10)%. (3.23) 

Of particular theoretical interest is the ratio of the branching ratios, defined as 

^C'Vn) . flel^. (3.24, 

in which a good part of the theoretical uncertainties cancels. This suggests that 
a future measurement of R{d'j/s'y) will have a large impact on the CKM phe- 
nomenology. 

Varying the CKM-Wolfenstein parameters p and f] in the range —0.1 < p < 0.4 
and 0.2 < rj < 0.46 and taking into account other parametric dependences stated 
above, the results (without electroweak corrections) are 

6.0 X 10"*^ < B{B Xd-f) < 2.6 x 10"^ , 
0.017 < R{d-i/s-i) < 0.074 . 

These quantities are expected to be measurable at the high-luminosity B facilities. 
3.4.3 B XJ+l- 

The inclusive B — > Xsl^l~ decay will also be accessible at the B factories. In 
comparison with the B — > Xgj decay, it presents a complementary and also 
more complex test of the SM since different contributions add to the decay rate 
(fig. 1^). Because of kinematic observables such as the invariant dilepton mass 
spectrum and the forward-backward asymmetry, it is particularly attractive. It 
is also dominated by perturbative contributions, if one eliminates cc resonances 
with the help of kinematic cuts. 

Using heavy quark expansion, nonperturbative corrections scaling with 1 / ml 



and 1 / can be calculated quite analogously to those in the decay B — > Xg'y p4 
However, there are also on-shell cc resonances, which one has to take into account. 
While in the decay B X^'j (on-shell photon) the intermediate ip background 
for example, namely B ipXg followed hj X'7, is suppressed and can be 
subtracted from the B Xg'y decay rate (see section lO] ), the cc resonances 



show up as large peaks in the dilepton invariant mass spectrum in the decay 
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Figure 9: One-loop contributions to the decay B — ^ Xgl^l 



B — > Xsl^l~ (off-shell photon). However, these resonances can be removed by 
appropriate k;inematic cuts in the invariant mass spectrum: In the 'perturbative 
window', namely 0.05 < s = (jrii+i- /mbY < 0.25, theoretical predictions for the 
invariant mass spectrum are dominated by the purely perturbative contributions, 
and theoretical precision comparable with the one reached in the decay B —>■ Xgj 
is in principle possible. 

The present status of the calculation of the perturbative contributions is the 
following: the effective Hamiltonian relevant to B ^ X^l^l' in the SM reads 

AC ^° 

H,jf{B ^ X^tl-) = ^ \, Y: aif^) 0,{ii) , (3.25) 



Compared with the decay B (see ( p.6| )), the effective Hamiltonian ( |3.2 

contains in this case two additional operators: 



09 = ^^{s^,PLb){m). 

(3.26) 

It turns out that the first large logarithm of the form log(mfe/M) (M = mw) 
arises already without gluons, because the operator O2 mixes into O9 at one loop 
(the pair cc in O2 can be closed to form a loop, and an off-shell photon producing a 
/ / pair can be radiated from a quark line). This possibility has no correspondence 
in the B Xg'j case within the SM. Consequently, the decay amplitude is ordered 
according to 

Gf log(mfe/M) (a,(m,) log(m,/M))" (LL), (3.27) 
Gf log(m,/M) «,(mfc) {a,{mb) log(m,/M))" (NLL), (3.28) 



which should be compared with ( |2.2| ) and { \i.'d\ ). To technically achieve the re- 
summation of these terms, it is convenient to redefine magnetic, chromomagnetic 
and lepton-pair operators C7, Og, Og, and Oiq and the corresponding coefficients 



as follows |62|, K3| 



Of - = = tI^C, {t = 7, 10). (3.29) 
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This redefinition enables one to proceed according to the three calculational steps 
presented in section |], when calculating the decay amplitude In partic- 



ular, the one-loop mixing of the operator O2 with the operator Og'^'^ appears 
formally at order g^, after the reshufHings in ( p.29| ). 



The QCD calculation up to NLL precision can be found in [62, However, 
the LL term in the series accidentially turns out to be small. In order to reach 
the same accurary as in the case of the NLL prediction for B Xg'-f one has to 
include the NNLL order contribution in the B —>■ XJ~^l~ calculation. 

Large parts of the NLL calculation in the decay B Xs'y, reviewed in section 



3^ , can be taken over and used in the NNLL calculation within the decay B — >■ 
Xsl^l~ . However, the complete NNLL enterprise - following the standard three 
steps in the formalism of effective theories (see section - is a formidable task: 

• Step 1: In [^] the complete Step 1 up to NNLL precision was presented. 
The authors did the two-loop matching for all the operators relevant to 5 — > 
Xsl^l~ (including a confirmation of the B — >■ X^'^ NLL matching results of 



32], 0, |3^). The inclusion of this NNLL contribution already removes 
the large matching scale {iiw) uncertainty of around 16% present in the NLL 
prediction oi B Xsl^l~ . As usual the partonic decay width is normalized 
by the semileptonic decay width in order to get rid of uncertainties due to 
the fifth power in mb- 

R^^^- is) - 1 dVih XJ^n 



One finds the following partial NNLL prediction ||60| : 

13{B —>■ Xsl~^l )cut: sG[0.05,0.25] = 

= B{B^X,eV) j^^^ ds[R[J;,,,{s)+^i/rnlR{s)+6,j^.R{s)] 
= 0.104 [(1.36 ± 0.18,ca/e ) + 0.06 - 0.02] 10"^ 

= {lA6±0.19scale)lO-^ (3.31) 

Si/^2R(s) and (5i/^2i?(s) are the nonperturbative contributions discussed in 
section p73| . The quoted error in ( p.31| ) reflects only the fib scale uncertainty. 



This purely perturbative uncertainty should get significantly reduced by 
contributions within Step 3 of the NNLL program, namely the two-loop 
QCD corrections to the matrix element of the four-quark operators 02- 
The error due to the uncertainties in the input parameters and to other 



contributions was not estimated in [^, at this intermediate stage of the 
NNLL calculation. 

Step 2: The most important NNLL contribution from the three-loop renor- 
malization group evolution of the Wilson coefficients from the matching 
scale fiw to the low scale fib, namely the three-loop anomalous dimensions 
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corresponding to the mixing of the four-quark operators Oi {i = 1...6) into 
the dipole operators O7 and Og, can be taken over from the NLL calculation 
in the decay B Xg'j [^. However, the analogous three-loop anomalous 



dimensions corresponding to the mixing of the four-quark operators into 



the operator Og is missing. In I^Oj an estimate was made which suggests 
that the numerical influence of these missing NNLL contributions to the 
branching ratio of -B ^ XJ^l~ is small. 

Step 3: Within the NLL B — > Xg'y calculation the two- loop matrix ele- 
ments of the four-quark operator O2 for an on-shell photon were calculated 
3l[] and quite recently confirmed in ||3^. This calculation was extended 



m 



to the case of an off-shell photon |6T|, which corresponds to a NNLL con- 
tribution relevant to the decay B — > Xgl^l" ■ The calculation includes also 
that part of the corresponding gluon bremsstrahlung which is needed to 
cancel infrared and collinear singularities of the virtual corrections. If one 
includes also this NNLL piece in the partonic NNLL prediction for the decay 



B l^l , one gets pT 

fO.25 
/0.05 



/■0.25 

/ dsRlLM = (1.25±0.08,,,,e) X 10-5 (3.32) 

J0.05 ^ 



Again the only error given corresponds to the uncertainty of the low scale 
/ifo. As expected the inclusion of the two-loop virtual corections to the four- 
quark operator O2 has reduced this scale ambiguity from ±13% down to 



±6.5%. The authors of [^] also analyse the error due to the uncertainty in 
the input parameter mc/mf, and find an uncertainty of ±7.6% within the 
partonic quantity. 

Within the Step 3 of the NNLL calculation, the renormalization group 
invariant two-loop matrix element of the operator Og is not calculated yet. 
Because this contribution includes no logarithms, the scale dependence of 
the NLL prediction is not sensitive to this NNLL contribution. 

One could think that within this perturbative window at low s G [0.05,0.25], 
one is only sensitive to C7 which would be redundant information, since we already 
know it from the decay B Xg^. However, as was explicitly shown in |63| , 
one is also sensitive to the new Wilson coefficients Cg and Cio and interference 
terms in the low s regime with s = rrii+i- /ml G [0.05, 0.25] (see fig. |10| where the 
various perturbative contributions to Rquark (with NLL precision) are plotted). 

Together with the decay B ^^7, the inclusive B — > XJ~^l~ decay will make 
precision flavour physics possible, if one can also measure the kinematic variables 



in the B Xgl^l decay precisely. As was first advocated in ||64 
• the invariant dilepton mass spectrum 

dV{B ^ Xsl^r ) / ds, (3.33) 
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Figure 10: Comparison of the different short- distance contributions to Rquark{s) 
(NLL precision), from . 



the forward-backward charge asymmetry 

A{s) = dcosO d^T{B ^ Xsl^r) / dsdcosO sgn(cos6 

and the decay rate of i? ^ ^sj, 



(3.34) 



(3.35) 



determine the magnitude and also the sign of the three Wilson coefficients Cj, Cg, 
and Cio, and allow for a model-independent analysis of rare B decays. For the 
measurements of these kinematic distributions, however, high statistics will be 
necessary. 



4 Indirect search for supersymmetry 

Today supersymmetric models are given priority in our search for new physics 
beyond the SM. This is primarily suggested by theoretical arguments related to 
the well-known hierarchy problem. The decay B — > Xg'j is sensitive to the mech- 
anism of supersymmetry breaking because in the limit of exact supersymmetry, 
the decay rate would be just zero: 

i3(5^X,7) = 0. (4.1) 

This follows from an argument first given by Ferrara and Remiddi in 1974 pH]. In 



that work the absence of the anomalous magnetic moment in a supersymmetric 



22 



abelian gauge theory was shown. The necessary mechanism of supersymmetry 
breaking, however, is unknown and leads to a prohferation of free parameters in 
the (unconstrained) minimal supersymmetric standard model (MSSM). 

There are two types of new contributions to flavour changing neutral currents 
in the MSSM: CKM-induced contributions, which are induced by a charged Higgs 
or a chargino, and generic new contributions, which are induced by flavour mixing 
in the squark-mass matrix. The structure of the MSSM does not explain the 
suppression of flavour changing neutral currents which is observed in experiments. 
This is the essence of the well-known supersymmetric flavour problem. 

In the framework of the MSSM there are at present three favoured concrete 
supersymmetric models. They solve the supersymmetric flavour problem by a 
specific mechanism through which the sector of supersymmetry breaking commu- 
nicates with the sector accessible to experiments: in the minimal supergravity 
model (mSUGRA) [|^, supergravity is the corresponding mediator; in the other 
two models this role is fulfilled by gauge interactions (GMSB) |^ and by anoma- 
lies (AMSB) [Q. Furthermore, there are other classes of models in which the 
flavour problem is solved by particular flavour symmetries pDI . 



Flavour violation thus originates from the interplay between the dynamics of 
flavour and the mechanism of supersymmetry breaking. The model-independent 
analysis of rare B and K decays therefore can contribute to the question of which 
mechanism ultimately breaks the supersymmetry and will thus yield important 
(indirect) information on the construction of supersymmetric extensions of the 
SM. In this context it is important to analyse the correlations between the dif- 
ferent information from rare B and K decays. 

As was already emphasized in the introduction, inclusive rare decays, as loop- 
induced processes, are particularly sensitive to new physics and theoretically 
clean. Neutral flavour transitions involving third-generation quarks, typically 
in the B system, do not pose yet serious threats to specific models. However, 
despite the relatively large experimental uncertainties, the rare decay B —>■ Xs'j 
has already carved out some regions in the space of free parameters of most of the 
models in the classes mentioned above (see and references therein). Once 

more precise data from the B factories are available, this decay will undoubtedly 
gain efficiency in selecting the viable regions of the parameter space in the var- 
ious classes of models; this may help in discriminating between the models by 
then proposed. In view of this, it is important to calculate the rate of this decay 
with theoretical uncertainties reduced as much as possible, and general enough 
for generic supersymmetric models. 

While in the SM, the rate for B —>■ is known up to NLL in QCD, the calcu- 
lation of this decay rate within supersymmetric models is still far from this level 
of sophistication. There are several contributions to the decay amplitude: be- 
sides the W t-quark and the H t-quark contributions, there are also the chargino, 
gluino and neutralino contributions. In most of the phenomenological analyses of 
the decay B these nonstandard contributions were not investigated with 

NLL precision as the SM contribution. However, as has already been pointed 
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Figure 11: SM, charged Higgs and chargino contribution at the matching scale. 

out, the step from the LL to the NLL precision is also necessary in order to check 
the validity of the perturbative approach in the model under consideration. It 
is possible that the restriction of the parameter space of nonstandard models, 
based on an LL analysis only, proves to be invalid after the NLL analysis is com- 
pleted. Moreover, it was already shown in specific new physics scenarios that 
bounds on the parameter space of nonstandard models are rather sensitive to 
NLL contributions (see below). 

Nevertheless, within supersymmetric models partial NLL results are available. 
The gluonic NLL two-loop matching contributions were recently presented ||72| . 



A complete NLL calculation of the B — > Xg'y branching ratio in the simplest 
extension of the SM, namely the Two-Higgs-Doublet Model (2HDM), is already 



available |]35| , |73| . In the 2HDM of Type II (which already represents a good 
approximation for gauge-mediated supersymmetric models with large tan/5 where 
the charged Higgs contribution dominates the chargino contribution), the B 

is only sensitive to two parameters of this model, the mass of the charged 
Higgs boson and tanp. Thus, the experimental data of the decay B — >• Xg'j allows 
for stringent bounds on these two parameters which are much more restrictive 
than the lower bound on the charged Higgs mass found in the direct search at 
LEP. One also finds that these indirect bounds are very sensitive to NLL QCD 
corrections and even to the two- loop electroweak contributions (see p5 . 



In a specific supersymmetric scenario is presented, where in particular 



the possibility of destructive interference of the chargino and the charged Higgs 
contribution is analysed. The analysis has been done under two assumptions. 
First it is assumed that the only source of flavour violation at the electroweak 
scale is that of the SM, encoded in the CKM matrix (minimal flavour violation). 
Therefore, the analysis applies to mSUGRA, GMSB and AMSB models (in which 
the same features are assumed at the messenger scale) only when the sources 
of flavour violation, generated radiatively between the supersymmetry breaking 
scale and the electroweak scale, can be neglected with respect to those induced 
by the CKM matrix. The second assumption is that there exists a specific mass 
hierarchy, in particular the heavy gluino limit. Indeed, the NLL calculation has 
been done in the limit 

fig ~ 0{mg,mg,mt^) > /iiy ~ 0{mw,mH+,mt,m^,mtJ. (4.2) 

The mass scale of the charginos (x), and of the lighter stop (^2) is the ordinary 
electroweak scale fiw, while the scale fig is characteristic of all other strongly 
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Figure 12: Upper bounds on the lighter chargino and stop masses from B — ^ Xg'y 
data in the scenario ( [4.2| ) if a hght charged Higgs mass is assumed; for tan [3 = 2 



(three lower curves) and 4 (three upper plots) the LL, NLL-running and NLL 



results (from the top to the bottom) are shown (see text), from WA 



interacting supersymmetric particles (squarks and gluinos) and is assumed to be 
of the order of 1 TeV. NLL QCD corrections have been calculated up to first 



order in /iv^/z^g including the important nondecoupling effects []7i |. 

At the electroweak scale ^y/i the new contributions do not induce any new 
operators in this scenario. Thus, the only step in the new NLL calculation beyond 
the one within the SM is Step 1, the matching calculation at the scale where 
we encounter the two new CKM-induced contributions of the charged Higgs and 
the chargino (see fig. |1T]): 

Cnll{,1^w) = Cnll^I^w) + C^LLif^w) + Cnll^I^w)- (4.3) 

One finds |]7i| that, in this specific supersymmetric scenario, bounds on the 
parameter space are rather sensitive to NLL contributions and they lead to a 
significant reduction of the stop-chargino mass region where the supersymmetric 
contribution has a large destructive interference with the charged-Higgs boson 



contribution. In fig. ^ the upper bounds on the lighter chargino and stop masses 
from B — >■ data in the scenario of ( [4.2|) are illustrated if a light charged Higgs 
mass of mH± = 100 GeV is assumed. The stop mixing is set to \9i\ < tt/IO which 
corresponds to the assumption of a mainly right-handed light stop. Moreover, 
< 500 GeV and all heavy masses are around 1 TeV. For tan/3 = 2 and 4 the 
results of the LL, 'NLL running' and NLL calculations are given. The result of 
neglecting the new NLL supersymmetric contributions to the Wilson coefficients 
is labelled as "NLL running" and illustrates the importance of the NLL chargino 
contribution |j7^ . 



Quite recently, this minimal flavour violation scenario was refined and ex- 
tended to the large tan(3 regime by the resummations of terms of the form 
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a^tan^^^fS |]73|, |7B]. The stability of the renormahzation group improved per- 
turbation theory was reassured for this specific scenario: the resummed NLL 
results in the large tan(3 regime show constraints similar to the LL results (see 
also 

For example, it is a well-known feature in the mSUGRA model, that depending 
on the sign oi Af ^ (where At denotes the stop mixing parameter) the chargino 
contribution can interfere constructively {At ■ /i > 0) or destructively {At ■ /x < 0) 
with the SM and the charged Higgs contribution. Therefore, the scenario Af^ > 
within this model requires very heavy superpartners in order to accommodate the 
B — >■ Xg'j data. But also the case AfH < is constrained in the large tanf3 regime 
where the chargino contribution is strongly enhanced (for details see JTS], |7^, [f7[]). 

However, all these NLL analyses are valid only in the heavy gluino regime. 
Thus, these calculations cannot be used in particular directions of the parameter 
space of the above listed models in which quantum effects induce a gluino contri- 
bution as large as the chargino or the SM contributions. Nor can it be used as a 
model-discriminator tool, able to constrain the potentially large sources of flavour 
violation typical of generic supersymmetric models. A complete NLL calculation 
should also include contributions where the gluon is replaced by the gluino. 

The flavour nondiagonal vertex gluino-quark-squark induced by the flavour 
violating scalar mass term and trilinear terms is particularly interesting. This is 
generically assumed to induce the dominant contribution to quark flavour transi- 
tions, as this vertex is weighted by the strong coupling constant Qs- Therefore, it 
is often taken as the only contribution to these transitions and in particular to the 
B — > Xgj decay, when attempting to obtain order-of-magnitude upper bounds 
on flavour violating terms in the scalar potential []78| , [79[] . Once the constraints 
coming from the experimental measurements are imposed, however, the gluino 
contribution is reduced to values such that the SM and the other supersymmetric 
contributions can no longer be neglected. Any LL and NLL calculation of the 
B —>■ rate in generic supersymmetric models, therefore, should then include 
all possible contributions. 

The gluino contribution, however, presents some peculiar features related to 



the implementation of the QCD corrections. In ref. I^O] this contribution to the 
decay B Xs'j has been investigated in great detail for supersymmetric models 
with generic soft terms. The gluino-induced contributions to the decay amplitude 
for B Xg'y are of the following form: 

«,(m,) (a,(m,) log(m,/M))" {LL), (4.4) 

al{mb){as{mb)\og{mh/M)r {NLL). (4.5) 
In the matching calculation, all factors as, regardless of their source, should be 



expressed in terms of the running with five flavours. In ||80| it is shown that the 
relevant operator basis of the SM effective Hamiltonian gets enlarged to contain 
magnetic and chromomagnetic operators with an extra factor of and weighted 
by a quark mass rrib or rric, and also magnetic and chromomagnetic operators of 
lower dimensionality where the (small) factor rrib is replaced by the gluino mass. 
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Figure 13: Matching of gluino-squark box on new scalar operators (left frame); 
mixing of new (scalar) operators at one-loop (a) in constrast to the vectorial 
operators of the SM (b) who mix at two-loop only (right frame). 



Furthermore, one finds that gluino-squark boxes induce new scalar and tensorial 
four-quark operators, which are shown to mix into the magnetic operators without 
gluons. On the other hand, the vectorial four-quark operators mix only with an 
additional gluon into magnetic ones (fig. 0). Thus, they will contribute at the 
NLL order only. However, from the numerical point of view the contributions 
of the vectorial operators (although NLL) are not necessarily suppressed w.r.t. 
the new four-quark contributions; this is due to the expectation that the flavour- 
violation parameters present in the Wilson coefficients of the new operators are 
expected to be much smaller (or much more stringently constrained) than the 
corresponding ones in the coefficients of the vectorial operators. This feature 
shows that a complete order calculation is important. 



In ref. |80| the effects of the LL QCD corrections on constraints on super- 
symmetric sources of flavour violation are analysed. To understand the sources 
of flavour violation that may be present in supersymmetric models in addition to 
those enclosed in the CKM matrix, one has to consider the contributions to the 
squark mass matrices 

= ( ""1^^ "^h"" ] + (4.6) 



Ff,LL + DfLL Ff^LR 

Ff,RL Ff,RR + ^f,RR 



(4.7) 



where / stands for up- or down-type squarks. In the super CKM basis where the 
quark mass matrices are diagonal and the squarks are rotated in parallel to their 
superpartners, the F terms from the superpotential and the D terms turn out to 
be diagonal 3x3 submatrices of the 6x6 mass matrices -Mj-- This is in general 
not true for the additional terms ([4.6| ), originating from the soft supersymmetric 
breaking potential. Because all neutral gaugino couplings are flavour diagonal in 
the super CKM basis, the gluino contributions to the decay b —>■ sj are induced 
by the off-diagonal elements of the soft terms 'm^LLi ''^},RRy ^},rl- 

It is convenient to select one possible source of flavour violation in the squark 
sector at a time and assume that all the remaining ones are vanishing. It should 
be stressed that one already excludes any kind of interference effects between 
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Figure 14: Gluino-induced branching ratio B{B Xg'y) as a function of x = 
m'^/m'^, obtained when the only source of flavour violation is Slr,23 (see text). 



different sources of flavour violation in this way. Following ref. all diagonal 



entries in ^.l^ ^'^d rr^ m\ are set to be equal and their common value 
is denoted by m?. The branching ratio can then be studied as a function of 

c _ {^\LL)ij jr _ ("^d,_Rij)*i ON 
OLL,ij — 2 ' ORR,ij — 2 ' 

hR,ij = ±2 (« ^ j)- (4- 9) 

rrig 

The remaining crucial parameter needed to determine the branching ratio is 
X = m|/m?, where rrig is the gluino mass. Figures n and 15 show the LL QCD 
corrections to the gluino contribution. 

In these figures the solid lines show the QCD corrected branching ratio, when 
only 5i/j,23 or Sll,23 are nonvanishing. The branching ratio is plotted as a func- 
tion of X, using mq = 500 GeV. The dotted lines show the range of variation of 
the branching ratio, when the renormalization scale /i varies in the interval 2.4- 
9.6 GeV. Numerically, the scale uncertainty in B{B Xs'y) is about ±25%. An 
extraction of bounds on the 6 quantities more precise than just an order of mag- 
nitude or less, would, therefore, require the inclusion of NLL QCD corrections. It 
should be noticed, however, that the inclusion of the LL QCD corrections has al- 
ready removed the large ambiguity on the value to be assigned to the factor as{fi) 
in the gluino-induced operators. Before adding QCD corrections, the scale in this 
factor can assume all values from 0{mb) to 0(mvi/): the difference between the 
value of B{B — > Xg'j) obtained when taking as{mb) and that obtained when tak- 
ing as{mw) is of the same order as the LL QCD corrections. The corresponding 
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Figure 15: Gluino-induced branching ratio B{B Xg'y) as a function of x = 
m?/m?, obtained when the only source of flavour violation is 5^^,23 (see text). 



values of B{B — > Xs'y) for the two extreme choices of are indicated in Figs. ^ 
and ^ by the dot-dashed lines (/x = mw) and the dashed lines (yU = nib). The 
choice /i = mw gives values for the non-QCD corrected B{B Xgj) relatively 
close to the band obtained when the LL QCD corrections are included, if only 
Sll,23 is nonvanishing. Finding a corresponding value of fi that minimizes the 
QCD corrections in the case studied in Fig. |n|, when only 5lb.,2z is different from 
zero, depends strongly on the value of x. In the context of the full LL result, it 
is important to stress that the explicit factor has to be evaluated - like the 
Wilson coefficients - at a scale /i = 0{mb). 

In spite of the large uncertainties still affecting the branching ratio B{B — > 
Xs7) at LL in QCD, it is possible to extract indications of the size that the 5- 
quantities may maximal acquire without inducing conflicts with the experimental 
measurements (see [PU|). 

Finally, it should be emphasized that a consistent analysis of the bounds on 
the sfermion mass matrix should also include interference effects between the 



various contributions. For this issue we refer to a quite recent paper |81|, where 



the interplay between the various sources of flavour violation and the interference 
effects of SM, gluino, chargino, neutralino and charged Higgs boson contributions 
is systematically analysed. New bounds on simple combinations of elements of 
the soft part of the squark mass matrices are found to be, in general, one order 
of magnitude weaker than the bound on the single off-diagonal element 5l_r,23, 
which was derived in previous work [^, ^ by neglecting any kind of interference 
effects. 
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5 Direct CP violation in 6/s transitions 



Detailed measurements of CP asymmetries in rare B decays will be possible in 
the near future. Theoretical predictions for the normalized CP asymmetries of 
the inclusive channels (see [091 B3l BH) within the SM lead to 



r(i?^x,/,7) + r(5^x,/,-7) 

acp{B ^ X,7) ^ 0.6%, acp{B ^ Baj) ^ -16% 
when the best-fit values for the CKM parameters p5| are used. 



(5.1) 
(5.2) 



The leading partonic contribution to the CP asymmetries is given by 

10-2 



acp{B Xs/dl) 



lap 



;i.l7 X Im [C2C;] - 9.51 X Im [C^C;] 



(5.3) 



+0.12 X Im [C2CI 



9.40 X Im 



0.013 c; 

p — ir] 



The large coefficient of the second term in ( ^.3|) has triggered an attractive 
scenario in which an enhanced chromomagnetic dipole contribution, Cg, induces 
a large direct CP violation in the decay B Xs'~f. Such a possible enhancement 
of the chromomagnetic contribution would lead to a natural explanation of the 
phenomenology of semileptonic B decays and also of charm production in B 
decays . 



An analysis for the leptonic counterparts is presented in |86]. The normalized 
CP asymmetries may also be calculated for exclusive decays; however, these re- 
sults are model-dependent. An example of such a calculation may be found in 



CLEO has already presented a measurement of the CP asymmetry in inclusive 
6 — > 57 decays, yielding ISH] 



acp{B X,7) = (-0.079 ± 0.108 ± 0.022) x (1.0 ± 0.030) 



(5.4) 



which indicates that very large effects are already excluded. 

Supersymmetric predictions for the CP asymmetries in B ^ Xg/dJ depend 
strongly on what is assumed for the supersymmetry-breaking sector and are, thus, 
a rather model-dependent issue. The minimal supergravity model cannot account 
for large CP asymmetries beyond 2% because of the constraints coming from 
the electron and neutron electric dipole moments |8^. However, more general 
models allow for larger asymmetries, of the order of 10% or even larger ||9^, 
|5B| . Recent studies of the B — > Xd'j rate asymmetry in specific models led to 
asymmetries between —40% and +40% or —45% and +21% In general. 



CP asymmetries may lead to clean evidence for new physics by a significant 
deviation from the SM prediction. From (^.2|), it is obvious that a large CP 
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asymmetry in the B Xg'j channel or a positive CP asymmetry in the inclusive 
B — Xd'j channel would be a clear signal for new physics. 

In [0] it was pointed out that the exclusive and inclusive decays of the form 



6 — > S7 and b dj, as well as their leptonic counterparts, provide a stringent test, 
if the CKM matrix is indeed the only source of CP violation. Using U-spin, which 
is the SU (2) subgroup of flavour SU{3) relating the s and the d quark and which is 
already a well-known tool in the context of nonleptonic decays |^4|, ^ , one derives 
relations between the CP asymmetries of the exclusive channels B^ ^ K*~'~f and 
B~ — ^ p~7 and of the inclusive channels B —* Xg'y and B Xd'y. One should 
make use of the U-spin symmetry only with respect to the strong interactions. 
Moreover, within exclusive final states, the vector mesons like the U-spin doublet 
{K*~ , p~) are favoured as final states because these have masses much larger 
than the (current-quark) masses of any of the light quarks. Thus one expects, 
for the ground-state vector mesons, the U-spin symmetry to be quite accurate in 
spite of the nondegeneracy of and nig. Defining the rate asymmetries (not 
the normalized CP asymmetries) by 

Ar(5- ^ V--i) = T{B~ I/-7) - T{B+ r+7) (5.5) 

one arrives at the following relation 



AT{B- ^ K*--f) + AT{B- ^ p-7) = be^c^exc (5.6) 

where the right-hand side is written as a product of a relative U-spin breaking 
bexc and a typical size Aexc of the CP violating rate difference. In order to give an 



estimate of the right-hand side, one can use the model result from fS^] for Aexc, 

Aexc = 2.5 X 10-' Tb. (5.7) 

The relative breaking bexc of U-spin can be estimated, e.g. from spectroscopy. 
This leads us to 

i(Mx* + mp) 

Certainly, other estimates are also possible, such as a comparison of fp and fx*- 
In this case one finds a very small U-spin breaking. Using the more conservative 
value for bexc, which is also compatible with sum rule calculations of form fac- 
tors (see [^), one arrives at the standard- model prediction for the difference of 
branching ratios 

\AB{B- K*--f) + AB{B- p--f) \ ~ 4 x 10"^ (5.9) 



Note that the right-hand side is model-dependent. Still (^.9|) is of some use, since 
a value significantly above this estimate would be a strong hint that non-CKM 
sources of CP violation are active. 

The issue is more attractive in the inclusive modes. Due to the l/rrib expansion 
for the inclusive process, the leading contribution is the free fe-quark decay. In 
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particular, there is no sensitivity to the spectator quark and thus one arrives at 
the following relation for the CP rate asymmetries p3| : 



AT{B ^ X,7) + AT{B ^ X,7) = 6,„eA,„,. (5.10) 

In this framework one relies on parton-hadron duality (besides in the long- 
distance contribution from up-quark loops, which is found to be rather small W^ )- 
So one can actually compute the breaking of U-spin by keeping a nonvanishing 
strange quark mass. However, it is a formidable task to do this for the CP 
asymmetries and it has not yet been done. The typical size of bmc can be roughly 
estimated to be of the order of |6j„c| ~ ml/ ml ~ 5 x 10^^; |Aj„c| is again the 
average of the moduli of the two CP rate asymmetries. These have been calculated 



(for vanishing strange quark mass), e.g. in and thus one arrives at 



\AB{B X,7) + AB{B Xdi)\ ~ 1 ■ 10"^ (5.11) 



Again, any measured value in significant deviation of (|5.11 ) would be an indica- 



tion of new sources of CP violation. Although only an estimate is given here, it 
should again be stressed that in the inclusive mode the right-hand side in ( p.ll| ) 



can be computed in a model-independent way with the help of the heavy mass 
expansion. 



6 Kl TT^uu and n^uu 

The rare decays n^uu and tt^uu represent complementary opportu- 

nities for precision flavour physics. They are flavour changing current processes 
induced at the one-loop level (see fig. 0) and are exceptionally clean processes. 
In particular, the —>■ tt^uu amplitude can be calculated with a theoretical 
uncertainty below 3% P7|] . 



This implies the important role of these decay modes for CKM phenomenol- 
ogy: they play a unique role among K decays, like the Bd ipKs mode among 
the B decays. They allow a measurement of one angle of the unitarity triangle 
without any hadronic uncertainties to a precision comparable to that obtained 



by the B^ — > i^Ks mode before the LHC era p8 



sm2/3=^ r^ = V^ -= (6.1) 



s 



where a is just related to the Wolfenstein parameter A = 0.22 via (1 — A^/2)~^; 
Pq{K~^) = 0.40 ± 0.06 is the internal charm contribution to tt^uu; this 

quantity is known up to next-to-leading QCD precision, and the dependence on Vtb 
is only of second order in A; Bi and B2 represent here the reduced branching ratios 
Bi = B{K+ 7r+z/z/)/(4.11 x lO^^^) and B2 = B{Kl 7rVz/)/(1.80 x lO^^"). 

The time-integrated CP violating asymmetry in B^ i'Ks is given by 
Acpi^i^Ks) = — sm2i3xd/ {I + x^) where Xd = Am/T gives the size of B^ — B^ 
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Figure 16: One-loop diagrams contributing to A' — > ttvv. 

mixing. With {sin2f3).,^yi, = {sin2(3)^Ksi obtains an interesting connection 
between rare K decays and B physics: 

^'■•(^■•«-^) - -Ac^m.)'-^, (6.2) 



which must be satisfied in the SM. As was stressed in all quantities in this 
'golden relation' (|6.2| ) - except for Pq(K^) - can be directly measured experimen- 
tally and the relation is almost independent of Vch- 

Besides their rich CKM phenomenology, the decays Kj^ — >■ -n^vv and 
-n^vv are very sensitive to new physics beyond the SM. In addition, the theoretical 
information is very clean and, thus, the measurement of these decays leads to very 
accurate constraints on any new physics model. Moreover, there is the possibility 
that these clean rare decay modes themselves lead to first evidence of new physics 
when the measured decay rates are not compatible with the SM. 

New physics contributions in Ki^ tt^uu and vr+i/z/ can be parametrized 

in a model-independent way by two parameters which quantify the violation of 
the golden relation ( |6.2| ) p9| , |100|| . New effects in supersymmetric models can get 



induced through new box diagram and penguin diagram contributions involving 
new particles such as charged Higgs or charginos and stops (fig. ^), replacing 



the W boson and the up- type quark of the SM (fig. p!q). 

In the constrained minimal supersymmetric standard model (MSSM), where 
all fiavour changing effects are induced by contributions proportional to the CKM 
mixing angles the golden relation ( |6.2|) is valid. Thus, the measurements of 
B{Kl — > vr^z/z/) and B{K^ — >• n^vv) still directly determine the angle /3, and 
a sigificant violation of ( p.2|) would rule out this model. 

At the present experimental status of supersymmetry, however, a model- 
independent analysis including also fiavour change through the squark mass ma- 
trices is more suitable. If the new sources of fiavour change get parametrized 
by the mass-insertion approximation, an expansion of the squark mass matrices 
around their diagonal, it turns out that SUSY contributions in this more gen- 
eral setting of the unconstrained MSSM allow for a significant violation of the 
golden rule. An enhancement of the branching ratios by an order of magnitude 
(in the case of —>■ vr+z/z/ by a factor 3) compared with the SM values is pos- 



sible, mostly due to the char gino- induced Z-penguin contribution |101|. Recent 



analyses [ |101| , |102| , |103[ within the uMSSM focused on the correlation of rare 
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Figure 17: Supersymmetric contributions io K nvv. 



decays and e'/e, and led to reasonable upper bounds for the branching ratios: 
B{Kl -K^vv) < 1.2 X 10"^°, and B{K+ n+uiy) < 1.7 x 10"^°. which should 
be compared with the latest numerical SM predictions ||104|| : B{K^ tt^uu) = 
(7.9 ± 3.1) X 10-^\ B{Kl vrVz/) = (2.8 ± 1.1) x IQ-^^ 

The rare decays n^uu and Kl tx^vv are specifically interesting 

in view of the suggested experiments at the Brookhaven laboratory (USA) ||108|| 
and at FERMILAB (USA) [|10|] [TTg and at KEK (Japan) [ml]. The current 
Brookhaven experiment E787 has already observed a single, but clean candidate 
event for K'^ n^uu in 1997 which corrsponds to the following branching ratio 
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B{K+ 



TX VV] 



;i.5i?:^) X 10- 



10 



(6.3) 



For the — > vr z/z/ mode, there is only an upper bound available ||106 



B{K, 



ttVz/) < 5.9 X 10 



-7 



(6.4) 



An indirect upper bound on B{Kl tt^uu), using the current limit on B{K~^ 
TT^vv) and isospin symmetry, can be placed |107| at 2 x 10~^. 



7 Summary 

In this paper we have reviewed the status of inclusive rare B decays, highlighting 
recent developments. These decays give special insight into the CKM matrix; 
moreover, as fiavour changing neutral current processes, they are loop-induced 
and therefore particularly sensitive to new physics. 

Decays modes such as -B — > and B —>■ Xsl^l~ (with specific kinematic 
cuts) represent laboratories for perturbative QCD. Nonperturbative contribu- 
tions play a subdominant role and they are under control thanks to the heavy 
mass expansion. The inclusive rare B decays are or will be accessible at the 
present e~ machines with their low background and their kinematic constraints 
(CLEO, BaBar, BELLE) and will make precision fiavour physics possible in the 
near future. 

Significant progress has been made during the last couple of years. The cal- 
culation of NLL (or even NNLL) QCD corrections to these decay modes has been 
performed. The theoretical uncertainty has been reduced below the 10% level. As 
was emphasized, the step from LL to NLL precision within the framework of the 
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renormalization group improved perturbation theory is not only a quantitative, 
but also a qualitative one, which tests the validity of the perturbative approach 
in the given problem. 

Inclusive rare B decays allow for an indirect search for new physics, a strat- 
egy complementary to the direct production of new (supersymmetric) particles, 
which is reserved for the planned hadronic machines such as the LHC at CERN. 
However, the indirect search at the B factories will imply significant restrictions 
for the parameter space of supersymmetric models and will thus lead to im- 
portant theoretically clean information for the direct search of supersymmetric 
particles. Within supersymmetric models the QCD calculation of the inclusive 
rare B decays has not reached the sophistication of the corresponding SM cal- 
culation. However, NLL analyses in specific scenarios already show that bounds 
on the parameter space of nonstandard models are rather sensitive to NLL QCD 
contributions. 

Detailed measurements of CP asymmetries in rare B decays will also be pos- 
sible in the near future. They will allow for a stringent and clean test if the CKM 
matrix is indeed the only source of CP violation. 

The rare kaon decays, n^uu and Kl —>■ n^uu, offer complementary 

opportunities for precision flavour physics. Besides the current Brookhaven ex- 
periment, several more are planned or suggested to explore these theoretically 
clean decay modes. 
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